Multi-wavelength anomalous diffraction at high x-ray intensity 
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The multi-wavelength anomalous diffraction (MAD) method is used to determine phase information in x-ray 
crystallography by employing dispersion corrections from heavy atoms on coherent x-ray scattering. X-ray 
free-electron lasers (FELs) show promise for revealing the structure of single molecules or nanocrystals within 
femtoseconds, but the phase problem remains largely unsolved. Due to the ultrabrightness of x-ray FEL, samples 
experience severe electronic radiation damage, especially to heavy atoms, which hinders direct implementation 
of the MAD method with x-ray FELs. We propose a generalized version of the MAD phasing method at high 
x-ray intensity. We demonstrate the existence of a Karle-Hendrickson-type equation for the MAD method 
in the high-intensity regime and calculate relevant coefficients with detailed electronic damage dynamics of 
heavy atoms. Our results show that the bleaching effect on the scattering strength of the heavy atoms can be 
advantageous to the phasing method. The present method offers a potential for ab initio structural determination 
in femtosecond x-ray nanocrystallography. 

PACS numbers: 87.53. -j, 61.46.Hk, 41.60.Cr, 32.90.+a 



Determination of the 3D structure of proteins and macro- 
molecules is crucial to understand their biological functions 
at the molecular level. X-ray crystallography has been widely 
used for structural determination [1], but it suffers from two 
bottlenecks: the phase problem and growing high-quality 
crystals. The phase problem [2, 3] is a fundamental obstacle in 
constructing an electronic density map from x-ray diffraction. 
Multi-wavelength anomalous diffraction (MAD) [4-6] with 
synchrotron radiation is one of the major achievements to ad- 
dress this issue. X-ray free-electron lasers (FELs) [7] promise 
to have a revolutionary impact on molecular imaging [8, 9], 
overcoming the crystal bottleneck. The unprecedented high 
x-ray fluence provides a sufficiently large number of photons 
to enable structure determination from diffraction measure- 
ments of streams of single molecules [9-11] and nanocrys- 
tals [12, 13]. However, due to an extremely high fluence that 
is ^100 times larger than the conventional damage limit [14], 
samples are subject to severe radiation damage [15]. The ul- 
trashort x-ray pulses generated by x-ray FELs enable us to 
carry out "diffraction-before-destruction" within femtosecond 
timescales to suppress nuclear motion [9]. Nonetheless, elec- 
tronic damage [9, 16, 17] during femtosecond x-ray pulses 
is unavoidable, leading us to consider "diffraction-during- 
ionization" [18]. This electronic radiation damage is par- 
ticularly challenging when addressing the phase problem by 
anomalous dispersion [19], because heavy atoms as anoma- 
lous scatterers will be more ionized than other atoms during 
intense x-ray pulses. Therefore, it has been speculated that 
MAD would not be an applicable route for phasing in the pres- 
ence of severe radiation damage [1, 19]. Here we propose a 
high-intensity version of the MAD phasing method based on 
a detailed description of the electronic response at the atomic 
level. In contrast to the speculation, our results show that 
MAD not only works, but also that the extensive electronic 
rearrangements at high x-ray intensity provide a new path to 
phasing. We will then demonstrate that this approach is appli- 



cable to the phase problem in femtosecond x-ray nanocrystal- 
lography [12, 13], which is one of the most prominent topics 
in x-ray FEL applications. 

X rays mainly ionize inner-shell electrons and subsequent 
relaxation (Auger decay and fluorescence) fills the inner-shell 
vacancy. Therefore, sequences of photoionization and relax- 
ation can strip off many electrons after absorbing several pho- 
tons [20]. For heavy atoms that have more than two sub- 
shells, a vacancy in deep inner-shells causes several relax- 
ation steps in the cascade through the subshells, resulting in 
further electron ejections [21]. To simulate the electronic 
damage dynamics, we use the XATOM toolkit [18, 22], where 
all rates and cross sections are calculated within the nonrel- 
ativistic Hartree-Fock-Slater method and multiphoton elec- 
tronic dynamics is described by sequential one-photon pro- 
cesses with all possible electronic configurations. Figure 1 
depicts the time evolution of populations for several charge 
states of an iron (Fe) atom where 27,783 coupled rate equa- 
tions were solved. The photon energy is 8 keV and the fluence 
is 5x10'^ photons//im". The pulse duration is 10 fs full- 
width-at-half-maximum (FWHM) with a Gaussian envelope. 
In these conditions, the neutral Fe is completely depleted and 
high charge states such as Fe^''^ are substantially produced 
by the end of the pulse. The pulse-weighted charge state av- 
eraged over time is about H-12, demonstrating the severe elec- 
tronic damage incurred during the x-ray pulse. 

This electronic damage affects not only the coherent scat- 
tering atomic form factor but also its dispersion correction. 
Near an inner-shell absorption edge, resonant elastic scatter- 
ing causes the atomic form factor to depend on the photon 
energy co. 



m,co)=f iQ)+/ (CO) + if" (CO), 



(1) 



where Q is the photon momentum transfer. The XATOM 
toolkit has been extended to compute the dispersion correc- 
tion, /' + //" [22]. In Fig. 2, one can see remarkable changes 
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FIG. 1. (Color online) Population dynamics for several selected 
charge states of Fe during an x-ray pulse. The thin dotted line indi- 
cates the Gaussian pulse envelope. See the text for parameters used. 
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FIG. 2. (Color online) Dispersion corrections of atomic form factors 
for selected configurations of several charge states of Fe. 



of the dispersion correction for different configurations and 
different charge states of Fe. Both /' and /" have a singular 
position at the A'-shell edge, which is shifted to a higher O) 
by ~1 keV as the charge state increases. The plotted curves 
in Fig. 2 correspond to the configurations of the ground state 
and the single-core-hole state (except for the neutral Fe) for 
given charge states. Since the MAD phasing method is based 
on the dispersion correction of heavy elements, it is inevitably 
required to take into account the electronic damage dynamics 
and accompanying changes of the dispersion correction under 
intense x-ray pulses. 

In the MAD phasing method, the Karle-Hendrickson equa- 
tion [23, 24] represents a set of equations of scattering cross 
sections at several different wavelengths (photon energies). 
The molecular scattering form factor is separated into normal 
and anomalous scattering terms and the phase information can 
be derived from their interferences. In this Letter, we show 
that a Karle-Hendrickson-type equation exists in the high- 
intensity regime with extensive electronic damage on anoma- 
lous scatterers. 

Let P be any protein (or any macromolecule) whose struc- 
ture we want to solve by coherent x-ray scattering. Let H 
indicate heavy atoms and Nh be the number of heavy atoms 
per macromolecule to be considered. Note that P excludes H. 
Our assumption is that only heavy atoms scatter anomalously 



and undergo damage dynamics during an x-ray pulse. It is jus- 
tified by the fact that the photon energy of interest is near the 
inner-shell ionization threshold of heavy atoms and the pho- 
toabsorption cross section a of the heavy atom is much higher 
than that of the light atom for a given range of o. For exam- 
ple, CTpe/cTc ~ 300 at 8 keV and there is almost no dispersion 
effect on carbon (C) near this photon energy (see Fig. Sl(b) 
in Ref. [25]). The scattering intensity (per unit solid angle) is 
evaluated by [25], 



dQ. 



^c(n) / dtgit)l^Pj{t) 

J — OO J 



(2) 



where j denotes a heavy atom index and / indicates a global 
configuration index. The global configuration for Nh heavy 
atoms is given by / = (/j ,/2, • • • ,Inh)- Here Ij indicates the 
electronic configuration of the j-th heavy atom, which is lo- 
cated at position Rj. This electronic configuration provides, 
among other things, information on the charge state of the 
atom. Pi{t) is the population of the I-th configuration at time t. 
It is assumed that the heavy atoms are ionized independently, 
so the population of / is given by a product of individual pop- 
ulations, Pi{t) = n^"jP/^. (f). ^ is the x-ray fluence and ^(f) 
is the normalized pulse envelope. Then the x-ray flux is given 
by ,^g{t), which is assumed to be spatially uniform through- 
out the sample. C(i2) is a coefficient given by the polarization 
of the x-ray pulse. 

In Eq. (2), Fp{Q) is the molecular form factor for the pro- 
tein (without any dispersion correction) and our purpose is to 
solve its amplitude and phase, Fp{Q) = |Fp (Q)| exp[/0^(Q)]. 
fij (Qi (o) is the atomic form factor (with the dispersion cor- 
rection) of the j-th heavy atom in its Ij-th configuration. It 
is most instructive to consider only one heavy atomic species. 
We introduce a molecular form factor for undamaged heavy 
atoms, 

Nh 

;=i 

where (Q) indicates the normal scattering atomic form fac- 
tor for the ground-state configuration of the neutral heavy 
atom. 

Now Eq. (2) can be expanded to demonstrate the existence 
of a generalized Karle-Hendrickson equation [25], 



dHQ,co) 



dQ. 



^C(n)[|FO(Q)| +|FO(Q)|'a(Q,to) 
+ |^^(Q)||^';j(Q)|MQ,«)cos((/.«(Q)-0«(Q)) 
+ I^^AQ)! \F^{Q)\ciQ,co)sm (./.^(Q) - ^{Q)) 
+ NH\f^{Q)f{a{Q,C0)-a{Q,(0)}], (4) 

where the new MAD coefficients depending on Q and (O are 



defined by 
a{Q,(o) 

c(Q,a)) 
a(Q,co) 



{m)y /« 

^p«K(Q)+4H} 



1 



(5a) 
(5b) 
(5c) 
(5d) 



Here Ih indicates the electronic configuration of the heavy 
atomspecies andP/^ = J"^c/fg(f)f/^(r) is the pulse-weighted 
averaged population for the ///-th configuration. The new 
MAD coefficients from Eq. (5a) to Eq. (5d) are atom-specific 
and must be calculated with electronic damage dynamics and 
configuration- specific atomic form factors. The coefficient a 
is an incoherent average of p with The coefficients b 
and c are the real and imaginary components of the averaged 
atomic form factor, respectively. The coefficient a in Eq. (5d) 
is obtained through a dynamical form factor defined by 



/H(Q,a),r) = £P/„(0//„(Q,a)), 



(6) 



which is a coherent average of the configuration-specific form 
factors over Ih at a given time t. This a coefficient thus rep- 
resents the effective scattering strength of the heavy atom. In 
contrast to the original Karle-Hendrickson equation, Eq. (4) 
is separated into light atoms (P) and heavy atoms (H) because 
both electronic damage and anomalous scattering are treated 
exclusively on H. If only the ground-state configuration is 
considered, i.e., no electronic damage occurs, then a = a and 
Eqs. (4) and (5) are reduced to the original Karle-Hendrickson 
equation except for the separation of P and H. 

This generalized Karle-Hendrickson equation constitutes a 
set of equations with different co at every Q. In Eq. (4) there 
are three unknowns: |F^(Q)|, |/^^(Q)|, and (j)°{Q} - 0° (Q) 
for a given Q. With three or more different (O, those un- 
knowns can be solved by the least-square method [26, 27]. 
Combined with Patterson or direct methods [28, 29], the am- 
plitude and phase of heavy atoms can be determined, so two 
unknowns of |/^^(Q)| and (j)p{Q) are to be solved with two 
different co. Once all amplitudes and phases of P and H are 
determined, it is straightforward to construct the total struc- 
ture of T = P + H. To obtain non-trivial solutions from the 
least-square method, the contrast between the coefficients at 
two different (O must be non-zero. This condition is fulfilled 
even in the presence of severe electronic damage as shown 
in the following discussion. We emphasize that, if the MAD 
coefficients are predetermined experimentally or theoretically, 
then one can solve the structure (amplitude and phase) from 
diffraction measurements directly, without any iterative phase 
retrieval algorithms [30]. 

Let us consider Fe atoms embedded in a protein and then 
radiate an x-ray pulse of 2x10'^ photons and 10 fs FWHM 




-0.05 
-0.10 
-0.15 
-0.20 
-0.25 
-0.30 
-0.35 



(c) c(0,fQ) Undamaged ■ 
A=12nm^ ■ 
A=4 ^m^ ■ 
^ s A=1 jim^ ■ 
A=Q,A iim^ ■ 



T 1 r 

(d) a(0,o))-a(0,(i!) 



Photon energy (keV) 



7 8 9 10 

Photon energy (keV) 



FIG. 3. (Color online) Coefficient.s in the generalized Karle- 
Hendrickson equation for Fe as a function of the photon energy. The 
fluence is given by 2x lO'- photons/A where A is the focal spot area. 



into the sample. Figure 3 displays a, b, c, and (a — a) for the 
forward direction (Q ~ 0) computed by the extended XATOM 
toolkit [22]. The fluence ^ is given by 2xl0'^ photons/A 
where A is the focal spot area. When high charge states are 
generated by ionization dynamics, the scattering strength is 
lowered due to the reduced number of scattering electrons and 
the change of the dispersion correction. The degree of lower- 
ing in a and b shows different behaviors below and above the 
neutral Fe edge. Below the edge, the scattering strength is 
less lowered than above the edge because ionization dynam- 
ics are dominantly initiated by L-shell ionization whose cross 
section is 8 times smaller than that of /iT-shell ionization [31]. 
Above the edge, /(T-shell ionization channels are open and lead 
to further cascade decays, stripping off more electrons. As a 
result, a and b are dramatically bleached out and their min- 
imum is deepened and broadened. The absolute value of c, 
which corresponds to the averaged absorption cross section, 
is decreased as the fluence increases. 

This bleaching effect on the scattering strength is benefi- 
cial to the phasing problem in two ways. First, the contrast 
of the coefficients to be exploited in the MAD method is en- 
hanced. Even though the scattering strength is lowered for 
all (0, Figs. 3(a)-(c) clearly show that the low intensity cases 
(long-dashed and short-dashed lines) display a contrast sim- 
ilar to the conventional MAD method (solid lines). For the 
high intensity cases (dotted and dash-dotted lines), the con- 
trast in d and b becomes even larger when (O is chosen be- 
low the edge and around the minimum. The contrast in c is 
reduced to some extent but is not completely eliminated. It 
is worthwhile to note that broadening of the edge at high in- 
tensity makes precision of (O less important in experiments. 
Second, it brings an alternative phasing method similar to 
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single isomorphic replacement (SIR) [3] or radiation-damage 
induced phasing (RIP) [32]. By choosing one (O below and 
one CO above the edge, one can create two datasets that differ 
only in the scattering strength of the heavy atoms, and then 
solve the rest of the structure by density modification. In this 
method, there is neither atomic replacement in sample prepa- 
ration like SIR nor chemical rearrangement during the x-ray 
pulses like RIR Therefore, the rest of the structure remains 
invariant in the two different datasets. 

Now we discuss experimental implementation of the gener- 
alized version of the MAD phasing method. We used a Gaus- 
sian pulse in the above calculations. However, when the scat- 
tering strength of dl/dQ. is measured at a particular Q and 
CO, the pulse shape g{t) and the fluence ^ at a given posi- 
tion X in the x-ray beam may differ from shot to shot. Fol- 
lowing the procedure in Ref. [33], we numerically confirmed 
that {dI[^{x)g{t)]/dD.) w dI[{^{x)g{t))]/dD. to within 3%, 
where ( • ) denotes an ensemble average. Then the total sig- 
nal can be obtained by integrating over the interaction vol- 
ume, / d^xdl[{,^{x)g{t))] /dQ. X «n,oi(x), where «moi(x) is 
the molecular number density. In this process, the MAD coef- 
ficients are calculated with given {JFg{t)), and the basic struc- 
ture of Eq. (4) remains unchanged. 

In our model, resonant absorption processes and shake- 
up/shake-off processes [34] are neglected. They tend to gen- 
erate further high charge states, so the contrast enhancement 
and the bleaching effect would remain after inclusion of these 
processes. We note that the effect of impact ionization [35] 
on coherent diffractive imaging may be suppressed by using a 
sufficiently short x-ray pulse [18]. 

The structure of the generalized Karle-Hendrickson equa- 
tion [Eq. (4)] can be fully functional for phasing of nanocrys- 
tals, which are of current interest for structural determination 
with x-ray FELs [12, 13]. In Eq. (3), contains the structure 
factors of the heavy atoms. In the case of crystals, the heavy 
atoms are regularly located and can contribute to the Bragg 
peaks when satisfying Q • (R, — Rj) = Inn (n: integer) for all 
/ and j. In Eq. (5d), a is expressed with the coherent aver- 
age over configurations, and ie^^i^ 
with the coherent summation over heavy atoms 



fore. 



f^l^ from Eq. (3) is expressed 
There- 

2 



P^dtgit) ///(Q,co,r)lJ^iexp[/Q 



R, 



plies that all heavy atoms are described by the same dynam- 
ical form factor This term is then responsible for the Bragg 
peaks (oc A^^). On the other hand, the term Nh \f^\^{a — a) rep- 
resents fluctuations from all different configurations induced 
by electronic damage dynamics, corresponding to the diffuse 
background Nh)- As shown in Fig. 3(d), (a — a) increases 
as the fluence increases. However, it is an order of magnitude 
smaller than a and not confined to the Bragg peaks, implying 
that the high x-ray intensity does not fully destroy the coher- 
ent signals. 

In conclusion, we have proposed the MAD phasing method 
in extreme conditions of ionizing x-ray radiations. We assume 
that the scattering factors of the light atoms of the protein do 
not vary significantly over the measured range of x-ray fre- 



quencies, and that these atoms have normal scattering and 
no ionization. It is also assumed that the heavy atoms are 
ionized independently and only one type of heavy atoms is 
considered. We believe that the method should work even 
if these assumptions are removed, because the most impor- 
tant consequence of high-intensity x-ray irradiation — multi- 
ple ionization of the heavy atomic species — has been fully 
taken into account. We have combined electronic response 
at the atomic level and molecular imaging during intense x- 
ray pulses, and demonstrated the existence of a generalized 
Karle-Hendrickson equation for the MAD method at high x- 
ray intensity. The relevant coefficients to be used in the MAD 
method have been formulated and calculated with damage dy- 
namics and accompanying changes of the dispersion correc- 
tion. We have shown that the generalized equation is still 
applicable to the phase problem even in the presence of se- 
vere radiation damage. The bleaching effect on the scatter- 
ing strength of heavy atoms, which unexpectedly enhances 
the coefficient contrast in the MAD method, can be beneficial 
to phasing. Our study opens up a new opportunity of solving 
the phase problem in femtosecond nanocrystallography with 
x-ray FELs. 
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